THEORY OF PASSIVE MAGNETIC FIELD TRANSPORT 



KRISTOF PETROVAY 

Eotvos University, Department of Astronomy 
Budapest, Ludovika ter 2, H-1083 Hungary 



Abstract. In recent years, our knowledge of photospheric magnetic fields went through a thorough 
transformation — nearly unnoticed by dynamo theorists. It is now practically certain that the 
overwhelming majority of the unsigned magnetic flux crossing the solar surface is in turbulent form 
(intranetwork and hidden fields). Furthermore, there are now observational indications (supported 
by theoretical arguments discussed in this paper) that the net polarity imbalance of the turbulent 
field may give a significant or even dominant contribution to the weak large-scale background 
magnetic fields outside unipolar network areas. This turbulent magnetic field consists of flux tubes 
with magnetic fluxes below 10 10 Wb (10 18 Mx). The motion of these thin tubes is dominated by 
the drag of the surrounding flows, so the transport of this component of the solar magnetic field 
must fully be determined by the kinematics of the turbulence (i.e. it is "passive"), and it can be 
described by a one-fluid model like mean-field theory (MFT). The recent advance in the direct and 
indirect observation of turbulent fields is therefore of great importance for MFT as these are the 
first-ever observations on the Sun of a field MFT may be applied to. However, in order to utilize 
the observations of turbulent fields and their large-scale patterns as a possible diagnostic of MFT 
dynamo models, the transport mechanisms linking the surface field to the dynamo layer must be 
thoroughly understood. 

This paper reviews the theory of passive magnetic field transport using mostly first (and occa- 
sionally higher) order smoothing formalism; the most important transport effects are however also 
independently derived using Lagrangian analysis for a simple two-component flow model. Solar 
applications of the theory are also presented. Among some other novel findings/propositions it is 
shown that the observed unsigned magnetic flux density in the photosphere requires a small-scale 
dynamo effect operating in the convective zone and it is proposed that the net polarity imbalance in 
turbulent (and, in particular, hidden) fields may give a major contribution to the weak large-scale 
background magnetic fields on the Sun. 
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1. Introduction 

1.1. "Passive" Fields vs. "Active" Fields: a Historical Review 

A basic rule of thumb in magnetohydrodynamics (MHD) tells us that the character of 
the interaction between motions and magnetic fields in a (high plasma beta) plasma 
is determined by the ratio of the Em magnetic and Ek kinetic energy densities. If 
Em <K Ek then the Lorentz force may be neglected in the equation of motion and 
our problem is reduced to the kinematical case. If, on the other hand, Em 3> Ek 
then the field will "channel" the flow and the only potential effect of the motions 
on the field is the generation of small- aplitude MHD waves: this is the strong field 
case. Finally, in the hydromagnetic case, when Em ~ Ek, there is a complicated 
interaction of flow and magnetic field. Of course we must be aware of the fact that 
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(B 2 ) > (B) 2 in general, so the Em total magnetic energy density may well exceed 
the energy density of the large-scale mean field. Besides, the validity of the above 
simple rule may possibly also be limited in two dimensions where a weaker magnetic 
field (consisting e.g. of a low filling factor set of strong sheets) could possibly also in- 
fluence the motion, owing to the topological constraint (the flow cannot "get around" 
the sheets). These latter points were recently brought into focus by Cattaneo and 
Vainshtein (1991). Nevertheless, apart from these rather obvious caveats, the simple 
rule summarized above can be considered as correct. This simple notion formed the 
background of MHD thinking in the 1950's and 60's when mean-field electrodynam- 
ics and mean-field MHD were developed (Parker, 1955, Steenbeck et at, 1966) for 
the treatment of the kinematic and hydromagnetic case, respectively. 

The picture however got more complicated in the period from the mid-sixties to 
the mid-seventies when solar observations (Sheelcy, 1966, Stenflo, 1973, Howard and 
Stcnflo, 1972) and numerical experiments (Weiss, 1964) showed that in the highly 
conductive turbulent solar plasma the magnetic field is concentrated into strong flux 
tubes with very little flux in between. The B t magnetic flux density inside the tubes 
is order of (or greater than) the B eq equipartition flux density defined by 

B 2 q /2p = pv 2 /2 (1) 

(SI formula; p is the permeability, p is the density, Vt = (u 2 ) 1 / 2 with v the turbulent 
velocity) . As a consequence of this realization, flux tube theory began to develop in 
the 1970's (see Parker, 1979, for a review of the results of this period). According to 
flux tube theory, the most important forces acting on a magnetic flux tube are the 
.Fd aerodynamic drag, the F m magnetic curvature force and the buoyancy; their 
approximate expressions are: 

2p Hp p R c a 

with Hp the pressure scale height, d the tube diameter and R c the curvature radius; 
in practice, one may put R c ~ I with I the characteristic scale of the turbulence. 
A comparison of these expressions shows that for a sufficiently thin flux tube, with 
a magnetic flux $ < $ cr = min {I 2 , Hp}B^ q /B^ , the drag will dominate and the 
surrounding flow will determine the motion, while thicker tubes may move more 
independently of the surrounding turbulence, under the action of dynamical forces. 
This implies that, for a given energy density (for a given / magnetic filling factor 
or D typical tube separation), the transport of the field may be either passive, i.e. 
fully determined by the flow between the tubes, or active, i.e. to a large degree 
independent of the flow in the non-magnetic component, depending on whether the 
field is organized in a large number of thin fibrils or in a small number of thick 
flux bundles. Our classification of magnetic field types from the point of view of 
transport should therefore be extended into a two-dimensional scheme: 
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(cf. also Fig. |l|). Traditional MFT, with its one-fluid approach, may be applied to 
passive fields, but it is clearly not applicable to the case of active fields as these 
would require a two-fluid description. (An interesting attempt at the construction 
of a more general, two-fluid MFT was made by Parker, 1982.) 



PW PM 




Fig. 1. Types of magnetic field structure from the point of view of field transport. Note that 
"passive" (P) does not imply "kinematical" : e.g. in case PM the presence of thin flux tubes in 
large numbers may significantly alter the flow pattern by posing a serious obstacle in the way of 
the flow. 

Now let us estimate the value of $ cr for the solar photosphere. Using Hp ~ 
100 km, I ~ 10 3 km (the granular scale) and assuming B t ~ B cq ~ 40 mT we find 
<& cr ~ 10 9 -10 10 Wb. (The assumption of B t ~ B cq is to some extent validated by 
the observations of Keller, 1994.) Flux tubes with such low flux values were rarely 
(if ever) observed on the Sun prior to the mid-1980's. As a consequence, all or most 
of the photosphcric magnetic field was thought to be in active form, and, from the 
late 1970's onwards, serious doubts arose concerning the relevance of MFT to the 
solar dynamo problem. The emphasis of research shifted to flux tube emergence 
calculations and numerical experiments. 

The observations of the last 5-10 years have brought a fundamental change in our 
perception of photosphcric fields; an important consequence of these new results may 
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be a partial "rehabilitation" of MFT as the basic theory of solar activity. Because of 
their importance, these recent observational results will be summarized separately 
in the next subsection. First however we should briefly deal with a number of other 
objections against MFT. 

Beginning with the famous paper by Parker (1975) on the "buoyancy problem", 
evidence has been mounting for a dynamo located at the bottom of the solar con- 
vective zone or below. It is then clear that the study of magnetic field transport is 
crucial in understanding the cause of this spatial restriction of the dynamo mecha- 
nism as well as in relating the fields seen at the surface to the dynamo deep below. 
In the past, MFT dynamo models usually ignored the transport effects, except tur- 
bulent diffusion (and occasionally a phenomcnological term imitating buoyancy). 
In fact, although some theoretical results concerning magnetic field transport were 
sporadically published (and will be cited in subsequent sections of this paper), this 
subject has only been reviewed once (Moffatt, 1983), and solar applications only 
began in earnest in the 1980's. 

MFT's usual limitation to the kinematical case with first order smoothing also 
used to be a subject of criticism. In recent years, however, a number of important 
papers on higher order and nonlinear effects (cf. Sects. 3.2-3.4) improved the position 
of MFT in this respect as well. 

In summary: although MFT is certainly no wonder-drug, it is a powerful and 
versatile theory which is able to reproduce some basic features of solar activity (e.g. 
the butterfly diagram) and which should be applicable to at least one important 
component of the solar magnetic field: to passive fields. 

1.2. The Impact of Observations 

Passive fields are freely advectcd by turbulence, so their fluctuating component, 
the turbulent field must have a characteristic scale or correlation length not exceed- 
ing that of the turbulence. Numerical experiments and turbulence closure models 
(De Young, 1980, Mencguzzi et al, 1981, Mencguzzi and Pouquet, 1989) actually 
show that in high Reynolds number non-helical turbulence the maximum of the 
magnetic energy spectrum lies at significantly (about an order of magnitude) higher 
wavenumbers than that of the kinetic energy spectrum. On this basis we may expect 
that the scale of the photosphcric turbulent field should be about 100 km, well below 
the resolution limit. How can we hope to observe this mixed-polarity field? 

The most straightforward approach is to try to increase the resolution of mag- 
netograph observations. The presently available resolution of ~ 1" is however only 
sufficient to give us a glimpse of the low-wavenumber end of the supposed turbulent 
magnetic energy spectrum. Observationally, pixels with a coherent motion lying far 
from the network and containing a small amount of flux above the noise level may be 
identified as the resolved component of the turbulent field (TF), called intranetwork 
(IN) field. Any component with polarities mixed on a scale smaller than the reso- 
lution element will remain invisible (hidden field, HF). After tentative detections in 
the 1970's (Livingston and Harvey, 1971), the systematic study of IN fields began 
in the 1980's (e.g. Zirin, 1986). The most thorough investigations to date were done 
by Martin (1990, 1994). These videomagnetogram studies indicate that the IN ele- 
ments move coherently from pixel to pixel, i.e. each of them is the manifestation of 
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a single flux tube (with a relatively large flux compared to the other, unseen tubes) 
rather than being a larger scale net polarity fluctuation in the distribution of smaller 
turbulent flux tubes. (This may however partly be a selection effect: the magneto- 
graph noise level determination method, as described e.g. by Murray, 1992a, does 
not distinguish between instrumental and physical noise, so the noise-dominated pix- 
els may actually also contain resolved net polarity fluctuations in the distribution 
of the smallest tubes.) In a typical quiet sun area, the unsigned flux density was 
found to be about 0.5 mT (Martin, 1990), comparable to the value for ephemeral 
active regions (ER) and for the mixed-polarity network (MN) together, as well as to 
the cycle- and latitude-averaged value for active regions (AR) and unipolar network. 
Noting that the small flux and high inclination of the smallest elements of the mixed 
network (Murray, 1992b), which contain most of its unsigned flux, indicates that 
this component of the solar magnetic field should also be classified as passive and 
its fluctuating component may be thought of as part of the turbulent field, from this 
alone we may conclude that turbulent fields must give a large fraction of the total 
unsigned flux density. 

But how much unsigned flux resides in the HF? Some indirect methods have 
been devised to study this problem. A rather firm la upper limit of 10 mT can 
be placed on (|B|) on the basis of the lack of any significant Zeeman broadening 
in unpolarized spectral lines (Stenflo and Lindegren, 1977). Transverse Zeeman 
linear polarization yields the difference in the field components perpendicular to 
the line of sight; again, only upper limits can be given, so the fields arc close to 
isotropy (Stenflo, 1987). Stenflo (1982) was able to place a lower limit of 1 mT on 
(|B|) on the basis of an analysis of linear polarization measurements, utilizing the 
Hanle effect. In a more detailed analysis of the Hanle effect Faurobert-Scholl (1993) 
found that the net unsigned flux density in the low photosphere is between 3 and 
6mT (a recent study, Faurobert-Scholl et at, 1994, yields a similar but somewhat 
lower value). This implies that the overwhelming majority of (|B|) is in the form of 
turbulent fields, in line with a theoretical conclusion based on vertical flux transport 
models (cf. Sect. 4.2 below). If the typical scale of the turbulent field is indeed not 
much below 100 km then the LEST telescope may actually enable observers to get 
a glimpse of a significant fraction of the total unsigned flux. 

With this large value of (|B|), a net surplus of one polarity as little as 1 or 2 
percent over scales exceeding a few times 10 4 km will yield a mean net magnetic flux 
density comparable to the large-scale background magnetic field measured outside 
unipolar network areas. Indeed, the only viable alternative to the possibility that this 
background field is due to the turbulent fields is to assume that the surplus magnetic 
flux in one polarity originates from the AR which decay into strong, active network 
elements and these elements never decay further into small, passive MN/IN/HF 
elements. Before the existence of IN fields and the smallness and high inclination 
of most MN elements were recognized, this assumption was indeed often quoted to 
explain the origin of the large-scale fields. New observations however indicate that 
TF are a more likely candidate to have the polarity surplus. A comparison of auto- 
and cross-correlations of the large-scale polar field by Stenflo (1992) has shown that 
the rotation law of individual elements is quite different from the rotation of the field 
pattern which must be constantly renewed on a timescale between 3 and 30 days; 
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the flux emergence rate in AR is not sufficient to replenish the mean field in such a 
short time. Furthermore, high sensitivity magnetograms of the solar disk now show 
that, in contrast to unipolar network areas, the line-of-sight magnetic field in weak 
field areas does not vary too much between center and limb (Zwaan, 1994), i.e. these 
fields have a considerable inclination that is difficult to understand if the polarity 
surplus resides in the strongly buoyant thicker network tubes. The weak large-scale 
field must therefore be the result of small net polarity imbalance in the low flux 
component of the MN and/or in the IN/HF, all of which are passive and cannot 
avoid turbulent reprocession into IN fields on a short timescale. (AR however may 
still be the ultimate source of this reprocessed flux.) Note that it is even possible 
that most of the large-scale background flux is due to the HF, the elements of which 
are presently unavailable to direct observations! 

2. The main transport mechanisms 
2.1. Basic Formalism 

Throughout this paper, we will apply the following notations. The letters i through 
q denote integers taking the values 1..D where D is the number of dimensions (2 
or 3). Other letters stand for reals. Boldface and overarrow denotes vectors, hat 
denotes second-order tensors. For repeated small integer italic indices summation is 
implied — but not for capital indices! d x = d/dx; d x — d/dx; di = d/dxf, Sij is the 
unit tensor, is the Levi-Civita tensor. 

We start from the first two Maxwell equations supplemented with Ohm's law: 

VxB = ^j VxE = -d t B (3) 

j = <T(E + uxB), (4) 
with t the time, u the velocity and a the (tensorial) conductivity. From this 

d t B = V x (u x B) - V x [(t)V) x B] , (5) 

or in components: 

dtBi = eijkdjCkimUiBk — eijkdjeum^ipd P B m . (6) 

In a gas, the r) = (^cj) -1 turbulent magnetic diffusivity is isotropic, so r\ki — V^ki 
and the induction equation finally takes the form 

d t ~B = Vx(uxB)-Vx ?j(V x B). (7) 

In a turbulent medium both u and B may be split into average and fluctuating 
parts: 

b=(B)+B' u = U + v. (8) 

(We take ensemble averages, ignoring here the problem of relating these expected 
values to actual measurements. The Reynolds rule (B') = (v) = is assumed.) Our 
aim is now to derive the evolution equation analogous to (0) for the (B) mean field. 
The average of (Q) is: 

t (B) = V x (U x (B) + £ ) - V x r?V x (B). (9) 
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The real task then consists in finding an expression for the 

£ = (v x B') (10) 

turbulent electromotive force. 

In a discussion of magnetic field transport, after the derivation of the equation for 
(B) we will face the problem of distinguishing the terms responsible for transport 
effects from other terms. To some degree this is a matter of definition: the only 
reason why we first wrote down the anisotropic form (^) of the induction equation 
is that (||) is the par excellence transport equation for an advected solenoidal vector 
field — its first term expressing advection with the flow, its second term describing an 
anisotropic diffusion. One is tempted to say that those terms of (|]) that are formally 
identical to the terms of (|J) should be referred to as "transport terms"; however, 
from a conceptual point of view, our choice of definition should be such that the 
appearance of transport effects and of non-transport effects, respectively, will not 
depend on the chosen reference frame — i.e. only full tensorial terms obeying tensor 
transformation rules should be referred to as "transport terms" . Our approach will 
be to define "transport effects" as described by those and only those tensorial terms 
of (^) (or of £) that in some reference frame(s) give rise to terms formally identical 
to those of (Q). 

In order to get an expression for £, we subtract (||) from ([?]) to get 

d t B' = V x (U x B' + v x (B) + G) - V x r?V x B', (11) 

G = v x B' £. (12) 

Now, assuming that the fluctuating magnetic (and velocity) fields have a correlation 
time r (independent of location) over which all memory is lost, we may write 

£ = {vx f d t B' dt) - (v x «9 t B')r. (13) 
Jo 

Substituting here (|ll|), we get an approximate expression for £. In principle, a more 
exact formula might be derived by leaving the integral in ( |l3| ) instead of using multi- 
plication of the integrand with r as a crude estimate. This would result in integrals 
involving Green's functions (Radler, 1976). However, for a practical application of 
the method an order-of-magnitude estimate must be introduced again at the end, as 
the correlations appearing in the integrands are not known. The Fourier transform 
method (Moffatt, 1983) encounters a similar problem (with the added difficulty of 
treating inhomogeneity in a Fourier representation) . It will spare us a lot of work 
to make this unavoidable approximation right here at the outset. 

After the substitution, a third-order correlation will appear in the expression for 
£ (itself a second-order correlation) owing to the quantity G. A similar expression 
for G, in turn, would involve a fourth-order correlation, etc., so an infinite hierarchy 
of moment equations appears — a problem familiar from turbulence theory. In order 
to make the problem tractable, this chain must be broken somewhere. This is 
usually achieved by the so-called smoothing, i.e. by assuming that all correlations 
higher than a certain order disappear. (An alternative method will be applied in 
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section 2.4). In the simplest case of first-order smoothing (also known as second 
order correlation approximation or SOCA) one assumes G = 0. A sufficient (but 
not necessary!) condition for this is that B' <C |(B)|, which is true if the Strouhal 
number v t r/l <C 1 and v t r/HB <C 1, Hb being the scale length of (B). (Another 
sufficient condition, of less importance in the solar context, is that the magnetic 
Reynolds number lv t /n <C 1.) The expression for £ will be derived in Sect. 2.3. 
First, however, we will discuss the general form of £ and the form of the transport 
equations. 

2.2. Symmetry Considerations and Double-Scale Analysis 

Contrary to most discussions of MFT, this paper concentrates on the transport 
effects and does not deal with the dynamo effects. (For a review of the MFT of 
the solar dynamo the reader is referred to the papers by Radler, 1980 and Hoyng, 
1994.) Dynamo terms may be defined as those tensorial terms that would lead to 
an increase of the field without limit (in the kinematic description), were there no 
other terms present. Unfortunately, transport terms and dynamo terms do not form 
two disjunct sets, so the two problems cannot be discussed independently in the 
general case. In order to avoid any disturbing dynamo effect in our "sterile" study 
of the transport, we assume that all dynamo terms vanish. A (strictly unproven but 
very likely) conjecture states that a sufficient condition for this is that the velocity 
field should be parity invariant in the sense that all its mean pseudoscalar quantities 
vanish. From this point on, we therefore restrict our discussion to the case of parity 
invariant velocity fields. 

For a discussion of what this implies for the flow, we introduce the term pre- 
ferred orientation. If a second order tensor with a non-degenerate eigenvalue can 
be constructed from the velocity field, then we say that the corresponding set of 
eigenvectors defines a preferred orientation. Furthermore, the flow will be said to 
have a preferred sense of motion, if a non-zero polar vector can be constructed from 
the velocity field. In the same spirit, a preferred sense of rotation is said to exist 
if a non-zero pseudo-vector functional f[u(x)] can be constructed and if the flow is 
either not parity invariant or f has the property = 0. A preferred direction im- 
plies either a preferred sense of motion or a preferred sense of rotation. A preferred 
direction also defines a preferred orientation — but the reverse does not hold. 

The u(x) velocity field must obey the equation of motion; from the structure of 
the equation of motion it follows that in the kinematical case any deviation from 
statistical isotropy must be the consequence of either a non-vanishing g gravity or 
of a non-vanishing uj large scale vorticity of the flow. The assumption of parity 
invariance then requires gw = 0. This shows that we must either limit our study to 
the equator, or ignore stratification, or ignore rotation. It is this last assumption, 
i.e. the neglect of large-scale rotation that seems to be the least restrictive in the 
solar case, as it still allows us to incorporate some indirect effects of rotation in our 
model, notably 

1. A deviation from spherical symmetry (but not from axial symmetry). In a 
spherical Sun with 9<pr polar coordinates we therefore allow dg ^ for mean 
quantities, but we require (8^) = 0. In this sense our discussion will be limited 
to 2-dimensional models (although the fluctuations may be 3D). As a particularly 
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simple case, we will also regard plane parallel geometry with the z-axis parallel 
to g and (d y ) — 0; this can be viewed as a local approximation to the spherical 
case (corresponding to r — > oo). 

2. The existence of other preferred senses of motion beside g. The parity invariance 
condition however implies that all such senses of motion must lie in one plane, 
which is why we had to assume axial symmetry. (Otherwise, with f, g, h being 
3 preferred polar vectors, the pseudoscalar (f x g)h would not vanish.) Any 
preferred sense of rotation must then be perpendicular to this plane. The plane 
will be the meridional plane. Then it follows that for any a polar tensor we have 
a^r = arc/) = a6ci> = &<j}9 = 0, while for pseudo-tensors agg = a rr = ag r = a r g = 0. 
(Otherwise, ag would define a polar vector outside the meridional plane or a 
pseudo-vector in the meridional plane, respectively.) 

3. In principle, a meridional circulation can also be included, as its large-scale 
vorticity would lie in the azimuthal direction where a preferred sense of rotation 
is not excluded. In the present discussion, however, for simplicity we assume 
U = 0, i.e. we neglect both rotation and circulation. 

Now if we assume that the scale Hb of the large-scale mean field is (much) larger 
than vtT then, within a correlation time t, mean quantities like £ will only "feel" 
the large-scale field within a radius small compared to Hb, so we can expand £ in 
terms of (B) as 

S i = -a ( £>{B k )-l3§>d i {B k ). (14) 

(The minus sign is just a convention.) From ([l3|) and ( |TT| ) it follows that each term 
of £ must include the Levi-Civita tensor as one of its factors. Therefore we may 
write 

P$=etii0$ ) + eatf$' ) +l%l ( 15 ) 

where and /3^ k will certainly vanish in SOCA. (This is because "there are not 
the proper number of indices" in the first moment of (|ll|). For the same reason, 
in SOCA we cannot expect higher order terms in the expansion Jli|).) We now 
split and into symmetric and antisymmetric parts, further splitting the 
symmetric part of cS T > into an isotropic and a "traceless" part: 

a\ k = ~aS lk + a lk + e ijk 7j (16) 

$P =Pik + eijkSj ( 17 ) 
an = &ij = &ji Pij = Pji (18) 

With this, the general expression of the turbulent electromotive force (neglecting 
the higher order fi\-' k term) is 

£ = a(B) - 5(B) — 7 x (B) - (J3V) x (B) - (8 x V) X (B), (19) 

or 

£i = aS lk (B k ) - a lk {B k ) - e ljk jj(B k ) - e ijk (3jidi(B k ) - S k d t (B k ). (20) 

(Note that the brackets are redundant and they only serve to emphasize that our 
definition of (3 and 5 differs from that used in some other papers.) 
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Writing this into (g) we see that the role of —7 is similar to that of U i.e. this 
term describes the advection of the mean field with an effective velocity —7 without 
real large-scale motions. This effect is called (normal) pumping of the magnetic field. 
7 being a polar vector, the pumping requires the existence of a preferred sense of 
motion of the field. Such a preferred sense of motion may e.g. be due to a gradient 
in p (density pumping), to a gradient in v t (turbulent pumping), to a topological 
asymmetry of the flow (topological pumping), etc. 

The role of the /3 term is analoguous to that of the diffusive term in The 
/3-term however does not only describe an anisotropic diffusion, but its form is even 
more general than the corresponding term in (^), including a "diffusion" of the 
magnetic field along its own direction. For this reason this term is sometimes split in 
two parts corresponding to the symmetric and antisymmetric parts of the diB k tensor 
(e.g. Radler, 1980); this is however just a formal manipulation as the coefficients of 
these terms will not be independent. 

In a parity invariant flow the pseudoscalar a will vanish. But what about a? 
Writing down the components of £ one finds that the role of this term shows some 
similarity to that of the 7-term: it apparently describes a pumping — with a sign that 
depends on the orientation of the magnetic field! If the direction of this anomalous 
pumping is characterized by a 7 vector then there exists a A vector perpendicular 
to 7 such that the field component parallel to A will be transported along 7, while 
the component perpendicular to A (and to 7) will be transported in the opposite 
direction. (In a similar way, the <5-term leads to an anomalous or "skewed" diffusion 
effect.) As all terms of £ must include the Levi-Civita tensor, the general form of a 
is 

otik = ^ijHjCik, (21) 

to be compared with the 7-term in (|l6|). The dependence on A must obviously be 
contained in the polar tensor C. Its general form is 

Cik = o~ik + bXiXk + ce klm \ m . (22) 

(A free factor before the first term was omitted, as it may be defined into 7.) As in 
the kinematical description the pumping may only depend on the orientation of the 
magnetic field, but not on its direction, C may only involve terms of even order in 
Xi, i.e. c = 0. The symmetry condition on C then implies b = —2. 

Summarizing these considerations, we find that the general SOCA expression of 
£ for a parity invariant flow must have the form 

£1 = -£ijil](o~ik - 2XiX k )(B k ) - e ijk jj(B k ) - £ i]k f3jidi(B k ) - S k di{B k ). (23) 

Substituting this into (^|) we find the general form of the mean field transport equa- 
tion in an axially symmetric geometry with Aj = 5i2, U = (and therefore 5 = 0): 

d t A^ = aeed^A^ + a dr dgd r A^ + a rr dlA^ + aedeA^ + a r d r A^ + aoA^ (24) 

(B e ) = -A+/r - d r A^ (B r ) = dgA^/r + cot QA^r (25) 
dt(B^) = beedl(B^)+b er dedr{B^)+b rr dl(B^)+b e de{B^)+b r d r (B^)+b {B^) (26) 
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where 

aee = bee = Pee/r 2 a er = bg r = 2f3 8r /r a rr = b rr = /3 rr (27) 

ae = (70 +le)/r+ {Pee cot 9 - Pe r )/r 2 (28) 

a r = (7r + 7r) + (ftr + + A?r COt 0)/f (29) 

«o = [(7r + 7r) + cot 0( 7fl + 70 )]/r - \p„ + Pee cot 2 + /3 9r (l + cot0)]/r 2 (30) 
&e = (7fl - 7e)/r + Pu cot 61 A 2 + drper/r + d e Pee/r 2 (31) 

6 r = (j r - %) + (P rr + &$$)l r + ®rPrr + dePer/r (32) 

fro = [(lr - It) + d g {ie - 7e) + drP^/r + d r (j r - ■%) (33) 
+ [cot 6 d e p^ - P u j sin 2 8)}/r 2 (34) 

The molecular diffusivity term was ignored, as it can always be defined into the 
/3-term and is usually negligible. Note that for non-parity-invariant turbulence the 
appearance of /3gj, and P^ r would give rise to dynamo terms, i.e. the corresponding 
part of the /3-term is a transport and a dynamo term at the same time, illustrating 
our above point on the impossibility of a clean distinction between dynamo and 
transport effects. (The dynamo effect in question was discovered by Radler, 1980, 
who mentioned it in the context of a "/3w-dynamo".) 

In the plane parallel case these equations reduce to the particularly simple form 

d t A y = p xx d 2 A y + 2p xz d x d z A y + p zz d 2 A y + { lx + %)d x A y + ( 7z + %)d z A y (35) 

(B x ) = -d z A y {B z ) = d x A y (36) 

d t {B y ) = p xx dl{B y ) + 2p xz d x d z {B y ) + p zz d 2 {B y ) + (d x p xx + d z p xz )d x {B y ) 

MdzPzz + d x p xz )d z (B y ) + [d z { lz - %) + d x { lx - %)]{B y ) (37) 

2.3. First Order Smoothing 

It remains to find the expressions of the coefficients appearing in the transport 
equations. For this purpose we substitute ( |IT| ) in ([l^) under the assumption U = 
G = and find 

£i = -Te ik i{VjV k )dj(Bi) - Te ik id J (v j v k )(Bi) 

+Te lkl (v k d J v j )(Bi) - Te jk i(v k d l v j )(Bi) (38) 

Comparing this with (^3|) and assuming that the fluid obeys the di{pu{) = anelastic 
continuity equation (a good approximation in the solar convective zone), taking 
A,: = SiQ and eiKQ / Owe arrive at 

Pij = T(viVj) Si = (39) 

j k = Td j (v j v k )/2 (40) 

7* - r^^dvKVidip/p + d K (v 2 K )/2) + (dj(vjv K ) - d Q (v Q v K ))/2 (41) 
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.1 _ (v K divi) . . 

kA * ~ • (42) 

We see that the normal pumping only contains a term proportional to the gradient 
of the turbulent velocity correlation tensor, i.e. it is a pure turbulent pumping. The 
anomalous pumping, in contrast, also contains terms related to the density gradient. 

In the transport equations ji and 7$ only occur in combinations 7$ ± 7*. In 
spherical coordinates these combinations are: 

70 + 70 = T^^{v0V r )d r p/p + T-^^(de(vj) + 2(v e v r ))/r 

+Td r (v e v r ) (43) 

7r + Ir = T^(vl)d r p/p + T^fcdrtf) 

+Td e (v e v r }/r + t((v 2 ) - (v$))/r (44) 

70 -70 = -^7^r(^0^)drP/p + T T:i | TT (9e(^) + 2(v e v r ))/r 

+r(^ r )/r + r cot 0((v e 2 ) - (v 2 ))/r (45) 



7r-7r= - T f^HFT< v ?) d rP// , + ' r i3T+T^-^) 

+r cot 0(« 0Ur )/r + r((w 2 } - ( w 2))/ r ( 46 ) 
In the plane parallel case we have 

7x+7* = T^^(v x v z )d z p/p + T^^d x {vl) + rd z {v x v z ) (47) 

7, +7, - r^{v 2 z )d zP /p + r^d z {v 2 z ) +rd x (v x v z ) (48) 
lx-lx = -T^^{v x v z )d z p/ 'p + T-^^d x (vl) (49) 

7, -% = -T^{vl)d z p/p + T-^d z (v 2 z ) (50) 

For horizontally isotropic turbulence the vertical component of the anomalous 
pumping clearly vanishes (as we can expect from symmetry considerations). For 
a horizontal anisotropy induced by rotational influence we expect (v 2 ) < {v 2 ,), 
1A2 ^ 2, 3 Al < 1 and (v r v e ) < 0; in this case, the density pumping (i.e. the 
first term in each of the above expressions, proportional to Vp) is directed down- 
wards and polewards for poloidal fields and in the reverse direction for toroidal 
fields. On the other hand, the anomalous part of the turbulent pumping (to which 
the main contribution is given by the second term) acts in the opposite direction, if 
the turbulent velocity decreases inwards (in the bulk of the solar convective zone). 
In reality, all the anisotropy parameters involved should be uniquely determined 
by the SI rotational velocity and by the form of the velocity covariance tensor for 
the "original turbulence" (i.e. for the non-rotating case). E.g. assuming a quasi- 
isotropic original turbulence, Kichatinov (1991) found that the density pumping is 
always perpendicular to the rotational axis. 
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2.4. Lagrangian Analysis for a Two-Component Flow 

Despite its generality and elegance, the first (and higher) order smoothing discus- 
sion of the magnetic field transport suffers from two major drawbacks. First, the 
conditions of its applicability are unclear (B' <C |(B)| and a small Strouhal number 
are just sufficient conditions); second, the transport effects found by this method 
do not readily lend themselves to a simple physical interpretation. An alternative 
method of finding an approximate expression for the turbulent electromotive force is 
the so-called Lagrangian analysis; in essence, this involves the neglect of the diffusive 
term in (|lT|), in contrast to SOCA, where the G-term was neglected. In the present 
paper we do not present a general, abstract Lagrangian analysis of the magnetic 
field transport problem (see Moffatt, 1983 for such a discussion in the incompress- 
ible case). Instead, we limit ourselves to the case of a very simple (but anelastic, 
not incompressible) plane parallel 2D flow. Because of its simplicity, our description 
will be able to give an easy-to-grasp physical picture of the transport. 

Our 2D flow is assumed to be periodic in the horizontal (x-) direction. Of each 
period of length 21, a fraction /„ is occupied by a local upflow, and a fraction fd is 
occupied by a downfiow (/„ + fd = 1). This flow pattern is persistent for a time r, 
then new up- and downflows will instantaneously form with a random phase shift 
compared to the old flow. At the same time, diffusion is assumed to smooth the 
large scale field, destroying any correlations with the flow. (See Moffatt, 1983, for 
the general background of these assumptions.) The vertical (z-) velocity is assumed 
to be horizontally constant in both the up- and downflows at a value of v u and v,i, 
respectively. The jump at the up/downflow boundary is of course unphysical, but 
the present model only serves to give us a clue to the essence of the transport effects. 
We introduce a large-scale horizontal magnetic field which may lie in the plane of 
the motions or perpendicular to it. The layer may be stratified in p, v = (v u — i>d)/2 
and (B). For simplicity, we only regard one type of inhomogeneity at a time. The 
horizontal velocity is found from the anelastic continuity equation for each of these 
cases. We further assume 

l~VT<^H (51) 

where H is the scale height of the inhomogeneous quantity on which the transport 
depends. Then we have (v%) <v 2 , so the kinetical energy flux may be written as 

2F kin /p = f u vl + fdvl (52) 

while the anelastic continuity equation implies 

fuVu + fdVd = 0. (53) 

In the case when there is no inhomogeneity in v, we expect F^m = 0, and the above 
relations imply f u = fd = 0.5. 

Before we discuss the individual transport effects, some general remarks have to 
be made. The first point is illustrated in Fig. ||. In a purely vertical flow, the vertical 
transport of the field does not depend on its angle to the horizontal: this shows that 
our results for a horizontal large-scale field will remain valid for a tilted field as long 
as no large-scale horizontal flows exist. A horizontal flow, on the other hand, will 
give a vertical lift to a tilted tube — this effect will be important in Sect. 2.6 below. 
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Fig. 2. The tilt of a flux tube does not influence its vertical transport in a purely vertical flow, 
but a horizontal flow will give it an extra lift. 

Another remark concerns the condition (|^). We will find that under this as- 
sumption the SOCA expressions for the field transport effects are returned. In the 
Sun H > I for any mean quantity (cf. Sect. 4), so this condition is certainly less 
restrictive than that of the small Strouhal number. Furthermore, the assumption 
B' -C |(B) | is not needed. In fact, even in SOCA it would have been sufficient to 
assume that it is only the part of the fluctuating field correlated with the turbulent 
velocity that is small, while large fluctuations of independent origin (e.g. due to tur- 
bulent concentration by small-scale turbulence) may exist. In other words, the thick 
lines in Figs. ||-^| may represent field lines of a large-scale diffuse field or thin flux 
tubes embedded in a non-magnetic plasma. In this latter case, the direction of the 
individual tubes is clearly also irrelevant for the transport, which implies that each 
of the transport terms appearing in the transport equation of (B) must also appear 
in the transport equation of (|B|), assuming that (|B|) is mostly due to a strongly 
intermittent, connected field structure, and not to a number of isolated small flux 
rings. The reverse of this statement is however not true, as we will see in Sect. 3.2 
below. 

2.5. Turbulent Diffusion 

For the study of turbulent diffusion we assume Vp = Vv = 0, while the large scale 
horizontal magnetic field may depend on z. With v x = and letting Bq(z) — 
B x (z;t = 0) we have B x (x, z;t) = Bq[z — v z {x)t\. Then, assuming vt <C Hb, we 
have 

S y = (v z B x ) = -Tv 2 d z B , (54) 

independently of the orientation of the field. The form of the coefficient, tv 2 is to be 
compared with the SOCA expression (|3^). The physical interpretation is of course 
straightforward: if e.g. the field increases downwards, upflows bring more flux than 
downflows, resulting in a net transport upwards. The use of a turbulent diffusivity 
or conductivity was first proposed by Sweet (1950) and Csada (1951). 

2.6. Density Pumping 

The only inhomogenous mean quantity is now p, increasing downwards. We study a 
single half-cell [—0.5 < x/l < 0.5], representing the whole pattern; x — is fixed at 
an up/downflow boundary with x > corresponding to a downflow. The horizontal 
flow is then 



Vx = 



(x + l/2)v/H p for x < 
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o : poloidal field lines perpendicular to the plane / scalar contaminant 
^ : flow velocity 



Fig. 3. Density pumping for a poloidal field. 



(l/2-x)v/H p for x > (55) 

Again, we assume vt <C H p . 

First we regard the case when the field is perpendicular to the plane of the 
motions, see Fig. |[ This corresponds to the case of "poloidal fields" if our 2D flow 
pattern is identified with the extreme case of rotation-induced anisotropy (meridional 
rolls or "banana cells"). The field then behaves like a passive scalar; for each mass 
element, freezing-in requires B/p =const. Neglecting the horizontal flow we have 
B(x, z; t)/ p{z) — Bo/p(z — v z t). (It is easy to see that the horizontal flow will only 
give a correction of order vt/H p .) With this 

£x = -(v z B y ) = ~B Tv 2 d z p/p (56) 

Writing this into (||) we see that the pumping is directed downwards. The interpre- 
tation is again simple: downmoving elements get compressed, and this amplifies the 
field in the downflows relative to the upflows. The role of the horizontal flow is to 
return the decreasing up- moving flux into the downflows. 




i=o t>o 



— : toroidal flux tubes 
: flow velocity 



Fig. 4. Density pumping for a toroidal field. 



Next we consider the case of a "toroidal" field. The horizontal flow will in this 
case lead to the appearance of a vertical component of the fluctuating field in a finite 
part of the flow, as shown in Fig. [|. For the angle one can work out 

tan9~4iJ p /2, (57) 
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to first order in vt/H p . With this 

A/2 

£ y = (v z B x - v x B z ) = Bo/l (v z + v x tan9) dx ~ -B rv 2 d z p/p. (58) 

J -1/2 

The pumping is directed upwards. The interpretation is the following. The constant 
vertical velocity implies that in the upflows and in the central parts of the downflows 
the flux density does not change. In the peripheral parts of the downflows, however, 
where the field is tilted, the horizontal flow gives the tilted field an extra vertical lift, 
resulting in a reduced flux transport in the downflow as compared to the upflow. 

The downwards directed pumping of a poloidal field by a 2D cellular flow was first 
pointed out by Drobyshevski (1977). Vainshtcin (1978) showed that for 3D isotropic 
turbulence the density pumping will vanish. Density pumping was recently studied 
by Kichatinov (1991) under the assumption of quasi-homogeneity of turbulence. 

t ill 



t=0 

: toroidal flux tubes 
: flow velocity 



Fig. 5. Turbulent pumping for a toroidal field. 



2.7. Turbulent Pumping 

The only inhomogeneous quantity is now v 2 . As discussed in Sect. 2.4, in this case 
an asymmetry between up- and downflows is to be expected. In the case of weak 
inhomogeneity, implied by (]5l]), we may use the linear relation 2.Fkin ~ — plvtdv 2 / dz . 
Comparing this with ( |52] ) we find Av — (v u + t>d)/2 ~ 2vAf ~ —Idv/dz and 
A/ = (/« — /d)/2 = Av/2v <C 1. Then all expressions can be expanded into a 
series to first order in A/; in this approximation i*kin/ P — 2t> 3 A/. 

In the case of a poloidal field no turbulent pumping is expected to first order in 
vt/H v 2 , as the effect of the horizontal flow is again of higher order, and we are then 
back to the case of the density pumping of poloidal fields. 

In the case of a toroidal field, with a manipulation similar to the previous cases 
(plus with an expansion in A/), one can work out the result: 

£ y = (v z B x - v x B z ) = -B rd z v 2 . (59) 

The pumping is directed downwards if the turbulent velocity decreases downwards. 
The interpretation of (one contribution to) the pumping is shown in Fig. 0: as a 
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consequence of the continuity equation and of a downward directed kinetical energy 
flux, the downfiows are faster and they process a larger amount of flux through the 
reference level than do the upflows, resulting in a net transport downwards. (Another 
contribution is due to the action of horizontal flow on tilted flux tube portions, an 
effect analoguous to the density pumping of toroidal fields.) 

Turbulent pumping was discovered by Zeldovich (1957), who christened the effect 
"turbulent diamagnetism" . This expression is sometimes still used, although it is 
in fact a misnomer. Real diamagnetism is a dynamical effect, resulting from the 
orienting effect of the magnetic field on the medium, while turbulent pumping is a 
purely kinematic phenomenon. There is not even a mathematical similarity between 
the two effects: diamagnetism exerts its influence through the \x permeability, i.e. it 
is formally more similar to turbulent diffusion than to the pumping that involves an 
advective-type term. 

Another confusion related to turbulent pumping is that its expression is often 
quoted in the form 

7 = V/3 ~ V(v 2 t). (60) 

In the case when Vr ^ this clearly differs from (|4^) and from ([59|). For 2D 
turbulence the first equality in ( |60| ) is indeed exact (Cattaneo et ah, 1988) and it is 
also correct in 3D in order of magnitude. However, we should remember that the 
expression (3 ~ tv 2 is itself an order-of- magnitude formula only, so its differentiation 
will result in huge errors. The second (approximate) equality in ( |60| ) should therefore 
be used with caution. 

3. Higher Order and Nonlinear Effects 

3.1. Geometrical and Topological Pumping 

Geometrical pumping relies on the non- vanishing higher order correlation (v x d x v z ). 
The appearance of such a correlation implies that either the upflows or the downfiows 
will have a mean horizontal divergence. The drag associated with these horizontal 
motions will then preferentially sweep the magnetic flux into the convergent compo- 
nent, which results in a net flux transport along the convergent component. For the 
sweeping effect to be significant in a correlation time r, a rather large (comparable 
to v z ) horizontal velocity component correlated with d x v z is needed, which requires 
a strong inhomogeneity by virtue of the continuity equation (Petrovay, 1993). Con- 
sequently, in turbulent flows this effect may only be important near a possible sharp 
boundary of the turbulent layer, and is therefore generally limited to low Reynolds 
number flows. 

Topological pumping relies on a topological asymmetry of the up/downflows: in 
a horizontal cross-section of the layer, upflow areas may be connected, with isolated 
downfiows or vice versa. A one-dimensional object like a field line will more easily be 
carried with the connected component. Of course, ultimately "what goes up must 
come down", but, given sufficient time, magnetic diffusivity may be able to smooth 
out the resulting field into a mean field concentrated towards the direction of the 
connected flow. Again, the role of diffusivity implies that the effect will be more 
important for low magnetic Reynolds numbers. 
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These effects were much studied and became undeservediy widely known in the 
seventies and eariy eighties (Drobyshevski et al., 1974, Drobyshevski et al., 1980, 
Arter, 1983a, 1983b, 1985) when it was hoped that they may offer a solution to 
the buoyancy problem pointed out by Parker (1975). The interest has dwindled 
after their unimportance was recognized in the high (magnetic) Reynolds number 
case. Recently, Petrovay (1991) used topological pumping as an argument against a 
supposed sharp lower boundary of the solar convective zone. 

3.2. Negative Diffusivity and Small-Scale Dynamo 

Nicklaus and Stix (1988) derived the expressions of a and $ for homogeneous turbu- 
lence in the fourth-order correlation approximation. One of their most interesting 
results is that the fourth-order correction to the diffusivity is negative, and it in- 
volves the helicity fluctuations. If the inhomogencity of the field is strong, i.e. if 
v t r -C Hb does not hold, then this negative contribution may overcompensate the 
normal diffusivity. 

A negative diffusivity implies that any small fluctuation will increase instead of 
being damped. In fact, no matter how large the overall scale of (B) is, it may 
always develop infinitesimal small-scale fluctuations. The dominating negative term 
in (3 for these small scales then means that these small fluctuations may reach large 
amplitudes before their growth is halted by nonlinear effects. This process does not 
influence the large-scale field, as the random fluctuations give no contribution to the 
net mean field; but it will influence (|B|). In effect, the consequence of the negative 
diffusivity at small scales will be the appearance of a source term in the transport 
equation for (|B|), beside the transport terms known from the transport equation 
of (B) (cf. Sect. 3.4 above). Turbulence closure computations (Dc Young, 1980, 
Leorat et al., 1981, Durney et al, 1993) as well as numerical simulations (Meneguzzi 
et al., 1981, Meneguzzi and Pouquet, 1989, Nordlund et al., 1992) indeed show this 
small-scale dynamo action. On their basis, in the high (magnetic) Reynolds number 
case the source term may be phenomenologically approximated by 

a t (|B|) = (|B|)[l-(|B|)/(/ m B cq )] Wt //, (61) 

where / m ~ 0.1 is the magnetic filling factor at the nonlinear saturation. 

Physically, the small-scale dynamo relies on the existence of small-scale helicity 
fluctuations in the flow which increase the unsigned flux density by stretch-twist-fold 
type motions. 

3.3. Magnetic Buoyancy 

Although buoyancy was much studied in the case of active fields, its description in 
MFT was for a long time limited to the use of phcnomcnological sink terms in the 
mean induction equation. Recently however a major breakthrough was achieved by 
Kichatinov and Pipin (1993). In their SOCA treatment of the problem they recover 
the expected result 7 oc (B) 2 for weak fields. For fields closer to equipartition 
however, the buoyant loss is significantly reduced by the Lorentz force. 



PASSIVE MAGNETIC FIELD TRANSPORT 



19 



3.4. Density Self-Pumping and "Turbulent Buoyancy" 

In Sect. 3 we have seen that anomalous pumping effects, e.g. density pumping, must 
rely on the anisotropy of the turbulence in a plane perpendicular to the pumping. 
The feedback of the magnetic field on the turbulence may itself induce such an 
anisotropy, thereby leading to the transport of the field. In particular, a horizon- 
tal anisotropy induced by a non-vertical magnetic field should lead to the vertical 
transport of the field. This problem was recently studied by Kichatinov and Rudiger 
(1992). For strong fields they found a pumping downwards, in line with the expec- 
tation (the field should encourage turbulence in rolls around the field lines, which 
corresponds to the "poloidal" case depicted in Fig. ||). Surprisingly, however, for 
weak fields the pumping was found to be directed upwards. The origin of this curi- 
ous phenomenon of "turbulent buoyancy", as called by the authors, is unclear. 

3.5. Other Effects 

The above list is far from being exhaustive, even if we only regard the transport 
effects studied to date. As we go to higher orders, the number of new effects increases 
(although their importance will decrease, unless the inhomogeneity is very strong). 
Other effects are mentioned e.g. in Krause (1978), Moffatt (1983), Nicklaus and Stix 
(1988). 

An interesting but neglected aspect of the transport problem is the role of vor- 
ticity in the magnetic field transport. The association of vorticity with narrow 
downflows and its interaction with flux tubes (Nordlund et al, 1992) suggests that 
a study of this problem would be worthwhile. 

4. Solar Applications 

4.1. General Assessment and Convection Kinematics 

Numerical simulations (Chan and Sofia, 1987) show that in a quasi-adiabatically 
stratified turbulent convective layer the / correlation length is close to the Hp pres- 
sure cale height. The scale heights of other mean quantities (p, v%, (B) etc.) are 
usually greater than Hp, so in the deep solar convective zone the condition I < H 
(though not I <C H) is well satisfied, and the effects of first order in l/H, discussed 
in Sect. 2, may be expected to dominate in the mean field equation (with the notable 
exception of the small-scale dynamo source term in the equation for (|B|), discussed 
in Sect. 3.2 above). On the other hand, near the boundaries of the convective zone 
(e.g. in the photosphere) we may have I 3> H and higher order effects may become 
important. 

But what is the relative importance of different first order effects? In particular, 
is there any region in the convective zone where the importance of dynamo terms is 
negligible, so that our discussion of "sterile" transport in a parity-invariant flow ap- 
plies? As the expressions of the terms appearing in £ involved velocity correlations 
of different orders, in order to answer this question a model of the solar convective 
zone (SCZ) with a reliable description of convection kinematics is needed. Unfortu- 
nately, the available models of solar convection were developed either for the needs 
of stellar evolutionists, who only require an average stratification (and perhaps the 
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extent of the overshoot) or for the needs of solar/stellar granulation observers, for 
whom nothing short of a full 3D numerical simulation is good enough. Dynamo the- 
orists, with their medium-level requirements are not well served by contemporary 
convection theory. Nevertheless, the last decade has brought a number of advances 
in the difficult problem of convection kinematics. This field would be worth a review 
in itself, so here we only recall some results more directly relevant to our present 
problem. 

The best available model of the SCZ that can be recommended for use by dy- 
namo theorists is the model by Unno, Kondo and Xiong (Unno et al, 1985, Unno and 
Kondo, 1989); in what follows: the UKX model. This is a fully non-local model; its 
most important assumption is that third-order correlations are proportional to the 
gradients of second-order correlations; this is tantamount to the assumption of weak 
inhomogeneity. More restrictive from the point of view of magnetic field transport 
theory is the assumption of a constant radial anisotropy ((vf.) = 2(vf) = 2(w?)). As 
we have seen in Sect. 2, the anisotropy of turbulence plays an important role in the 
transport. Petrovay (1992) computed the anisotropy of low Prandtl number turbu- 
lent convection in the homogeneous (Boussinesq) case. Inhomogeneity may however 
significantly modify the anisotropy. In a later work, Petrovay (1993) evaluated the 
full anisotropy under the (invalid) assumption that inhomogeneity only influences 
the part of horizontal velocity that is correlated with the radial motions. In a 
complementary work, Kichatinov and Riidiger (1993) computed the full anisotropy 
caused by inhomogeneity in the case of quasi-isotropic original turbulence; this lat- 
ter work also incorporates the effect of rotation. Although a complete solution of 
the anisotropy problem is still not around, the available results suggest that the 
anisotropy (both radial and horizontal) in the SCZ is significant, but not extreme. 
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Log depth (m) 

Fig. 6. Characteristic velocites as functions of depth in the solar convective zone. Note that only 
an the upper limit on a compatible with the given kinetic energy and an upper limit on the density 
pumping velocity corresponding to the unlikely case of maximal anisotropy (meridional rolls) are 
displayed. 

Fig. fi shows the characteristic velocities corresponding to different terms in the 
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mean field transport equation as functions of depth in the SCZ. The relative mag- 
nitudes of the velocities characterize the relative importance of these terms. It is 
apparent that with the exception of the bottom of the SCZ the a term is by orders 
of magnitude smaller than the transport terms. The characteristic velocity differ- 
ence due to the shear by differential rotation (not shown) has a similar behaviour. 
This implies that with the exception of the bottom of the SCZ the dynamo terms 
are indeed negligible in comparison with the transport terms: in a simplistic treat- 
ment the solar dynamo may be split into a thin generation layer at bottom and a 
deep turbulent zone above; the distribution of the field generated in the dynamo layer 
will be determined by a dynamical equilibrium of turbulent transport processes in the 
convective zone. Of course we should keep in mind that part of the magnetic field 
intruding into the convective zone from the dynamo layer below will be in active 
form (e.g. the thick flux loops giving rise to AR). The decay of such thick active 
tubes may act as a local source of passive fields in the convective zone. The am- 
plitude of this active source may be determined from observations of AR decay at 
the surface; but in deeper layers it presents a large uncertainty in our discussions of 
passive field transport. 

As mentioned above, higher order effects may come into play near the boundaries 
of the SCZ. This is likely to be the case in the photosphere where I 3> Hp, but 
unfortunately very little is known about flux transport in these circumstances. In 
an ongoing effort (Faurobert-Scholl et at, 1994) the observed photospheric height- 
dependence of (|B|) is compared to vertical transport equilibrium models (as yet 
only including first order effects). In agreement with observations, these models 
indicate that (|B|) lmT at the top of the photosphere. Whether such a low 
magnetic energy density is sufficient to explain the observed basal chromospheric 
CallK flux, as proposed by Durney et al. (1993), remains to be seen. 

The lower boundary of the SCZ presents an even more vexing problem. Phe- 
nomenological non-local mixing-length models traditionally predict a very sharp 
lower boundary with v% suddenly dropping to zero in a thin layer. As pointed out 
by Petrovay (1991), in this case the extremely low velocity scale height should lead 
to strong geometrical and topological pumping. The flow topology in this layer is 
expected to be rf-type in the notations of Petrovay (1990), so all flux would be re- 
moved upwards from the dynamo layer. This contradiction seems to make it unlikely 
that the SCZ has a sharp lower boundary. In the models, the sharp bottom is a 
consequence of the a priori assumption of a good correlation between velocity and 
temperature fluctuations which is not granted in the overshooting layer. Indeed, the 
more realistic UKX model shows a much more gradual decrease of at the lower 
boundary with H v 2 ~ Hp. This result agrees with numerical experiments (Ludwig, 
1994, Brandenburg, private comm.) and helioseismological constraints (Monteiro 
et al., 1994). The problem is further complicated by the influence of magnetic fields 
and differential rotation on the structure of this important layer of the Sun. 

4.2. Vertical transport 

Fig. |^ shows that the most important vertical transport effects are turbulent diffu- 
sion and turbulent pumping. The vertical distribution of passive magnetic flux in 
the SCZ should therefore be determined by a dynamical equilibrium of these pro- 
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cesses, possibly influenced by density pumping and/or an active source. As ft, has 
a maximum at a shallow depth of 200 km in the UKX model, we may expect that 
the downwards directed turbulent pumping will "press" the field down to the bot- 
tom of the convective zone as far as turbulent diffusion, trying to "smooth out" the 
distribution, will allow. If for simplicity we regard a ID plane parallel model with 
(By) = 0, the only non-trivial component of the transport equation (|35|)-(|36|) is 

= d z (/3 zz d z (B x )) + d z [( lz +%)B x }, (62) 

with coefficients given by (|39|)-(|4l|). 

The vertical passive magnetic flux transport through the SCZ in a formulation 
similar to (|6^) was apparently first considered by Drobyshevski (1977). By analogy 
with the density pumping of a poloidal magnetic field in 2D convection he proposed 
that the vertical transport is determined by a strong density pumping downwards. 
Soon thereafter, Vainshtein (1978) showed that vertical density pumping will vanish 
for horizontally isotropic turbulence, and Schiissler (1984) stressed the importance of 
turbulent pumping in the vertical transport. The confusion concerning the correct 
expression of turbulent pumping, discussed above (end of Sect. 2.7), however led 
many to the false belief that turbulent pumping is negligible in the SCZ. As a 
consequence of this, density pumping made a comeback in the first quantitative 
study of the problem by Spruit et al. (1987). The heuristic transport equation used 
in this work was in effect equivalent to (|6^) containing only a (downwards) density 
pumping term and turbulent diffusion. An approximate analytic solution was also 
presented. Kichatinov's (1991) study of the density pumping was applied to the SCZ 
by Krivodubskij and Kichatinov (1991). 

The somewhat speculative early models came into a different light with the new 
observational results summarized in Sect. 1.2 above. In order to see what these 
empirical constraints imply for the vertical distribution of passive magnetic flux, 
Petrovay and Szakaly (1993) solved equation ( |62| ) using boundary conditions taken 
from the observations and including turbulent pumping and turbulent diffusion, 
but neglecting the contribution due to the third and fourth terms in ([58]) . Petrovay 
(1994) extended this work by including the neglected terms and the effects of spheric- 
ity (in an approximation that still preserved the mathemtical one-dimensionality of 
the problem). Equation (|6^) can be integrated to yield 

Pzzd z (B x ) + ( lz +%)(B x ) = -L m , (63) 

where the L m integration constant is the net flux loss from the Sun and it must 
obviously lie between and Vt tS (B x ) s ('s' index refers to values near the surface, 'b' 
to those at the bottom of the SCZ). Observations show that outside unipolar network 
areas (B z ) is typically order of 0.01-0.1 mT; this may also be applied to B x at some 
depth below the surface, assuming that averaging over the solar surface the ratio 
\(B X )\/\(B Z )\ ~ 1. Figure ^| (left) shows the numerical solutions of the transport 
equation with different boundary conditions. It is apparent that the flux density 
increases by several orders of magnitude to the bottom of the SCZ, confirming our 
expectation. Turbulent pumping therefore has the important role of "helping" the 
dynamo by sweeping the escaping flux back to the dynamo layer very effectively. It 
is also apparent that a relatively low flux density of ~ 100 mT at the bottom will 
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lead to a surface field comparable to that observed. Indeed, it is quite possible that 
all the net flux observed on the surface originates from passive fields in the dynamo 
layer by direct diffusion. On the other hand, we must remember that active source 
terms (neglected in (|62|)) may also be a significant source of net field. The observed 
rate of flux emergence in AR is too low to significantly modify the distribution 
(dynamical equilibrium is reached on a timescale ~ 20 days, comparable to the 
turbulent correlation time in the deep SCZ!) but AR may still be an important source 
of flux on a long term. In this latter case, the flux input by AR would be reprocessed 
by turbulent transport through the whole convective zone and would reappear again 
and again in a passive form. Finally, the possibility cannot be excluded that a 
much stronger active source operates in the deep SCZ in the form of a large number 
of "failed active regions" (Petrovay and Szakaly, 1993): rising flux loops decaying 
before they could make it to the surface. 




Fig. 7. Mean net flux density (left) and mean unsigned flux density (right) as functions of 
depth in the SCZ. Short dashes denote B cq (left) and f m B eq (right); other curves are solutions 
with different boundary conditions (see Petrovay and Szakaly, 1993, for a detailed legend). 



An equation similar to (p2[) , with the added source term (|6l| ) should determine 
the equilibrium distribution of the unsigned flux density (\B X \). The integration of 
this equation (Petrovay and Szakaly, 1993, Petrovay, 1994) yields the results shown 
in Fig. (right). Note the weak dependence on the boundary conditions applied 
and note that the surface value of the unsigned flux density (2.5-7.5 mT) is not a 
boundary condition here but an output in nice agreement with the observations of 
Faurobert-Scholl (1993). In fact the results show that even the resolved unsigned 
passive flux density of 0.5-1 mT (cf. Sect. 1.2) is incompatible with a transport 
equation without source term; this can be regarded as indirect evidence for the 
operation of a small-scale dynamo mechanism in the Sun. 
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4.3. Horizontal and 2D Transport 

The idea that horizontal transport of passive fields by meridional circulation and/or 
turbulent diffusion may have an important role in the solar dynamo goes back to 
Babcock (1961) and Leighton (1964) who proposed that the poleward drift of /- 
polarity background fields at high latitudes is due to these effects. These two effects 
were combined into a single model by Sheeley et al. (1985) and further developed 
until a detailed comparison with synoptic magnetic maps yielded very positive results 
(Wang et al., 1989). Despite their success in interpreting the synoptic maps, these 
ID horizontal transport models seem to be in conflict with the results of Stcnflo 
(1992), discussed in Sect. 1.2, which imply that the large-scale field pattern should 
be refreshed on a timescale of 3-30 days. Indeed, as we have seen, the vertical 
transport models yield a vertical transport timescale of ~ 20 days: over this time, 
the flux is vertically processed through the whole of the SCZ. The contradiction may 
disappear if we assume that the horizontal transport processes regarded by Wang 
et al. (1989) act throughout the whole SCZ. This would give a natural explanation 
for the remarkably low diffusivity (600 km 2 /s, much lower than tv\ ~ 10 4 km 2 /s for 
the supergranulation) required by the comparison with the synoptic maps: this is 
about the value of tv% in the deep SCZ. 

In order to check whether the results of the horizontal diffusion models remain 
valid, detailed 2D (horizontal-!- vertical) models of magnetic flux transport are clearly 
needed. A first step in this direction was made by Wang et al. (1991) where the 
vertical dimension was included by a kind of vertical truncation of the equations. In 
Budapest we are presently working on a fully 2D model of the field transport; this 
model involves the numerical solution of equations (|24|-[2(]) with proper boundary 
conditions. Preliminary results from a fast implicit integration of a simplified, ap- 
proximate form of the poloidal equation on a uniform 200 x 1000 grid show that the 
behaviour of the solution qualitatively agrees with the expectation, i.e. turbulent 
pumping presses down the field. 

A potential use of 2D transport models is to forge a link between the field distri- 
bution seen on the surface (as represented by synoptic maps or spherical harmonic 
expansions) and the field in the dynamo layer, thereby imposing an observational 
constraint on dynamo models. But the transport may also turn out to have a more 
fundamental role than this. If the ID results of Wang et al. (1989) are confirmed 
in 2D, then the "poleward branch" of the butterfly diagram is fully due to field 
transport effects, instead of being a dynamo wave. In fact, Kichatinov (1991) has 
proposed that transport effects may interpret both branches of the butterfly diagram: 
the horizontal density pumping should transport the toroidal field equatorwards, the 
poloidal field polewards. Clearly, the full role of flux transport mechanisms in the 
dynamo can only be judged after full 2D transport models become available. 
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